Optimal estimation of entanglement in optical qubit systems 
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We address the experimental determination of entanglement for systems made of a pair of polar- 
ization qubits. We exploit quantum estimation theory to derive optimal estimators, which are then 
implemented to achieve ultimate bound to precision. In particular, we present a set of experiments 
aimed at measuring the amount of entanglement for states belonging to different families of pure 
and mixed two-qubit two-photon states. Our scheme is based on visibility measurements of quan- 
tum correlations and achieves the ultimate precision allowed by quantum mechanics in the limit of 
Poissonian distribution of coincidence counts. Although optimal estimation of entanglement does 
not require the full tomography of the states we have also performed state reconstruction using two 
different sets of tomographic projectors and explicitly shown that they provide a less precise deter- 
mination of entanglement. The use of optimal estimators also allows us to compare and statistically 
assess the different noise models used to describe decoherence effects occuring in the generation of 
entanglement. 
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I. INTRODUCTION 



The sort of quantum correlations captured by the no- 
tion of entanglement represents a central resource for 
quantum information processing. Therefore, the precise 
characterization of entangled states is a crucial issue for 
the development of quantum technologies. In fact, quan- 
tification and detection of entanglement have been exten- 
sively investigated, see [ll-Q for a review, and different 
approaches have been developed to extract the amount of 
entanglement of a state from a given set of measurement 
results [3-0|- Of course, in order to evaluate the entangle- 
ment of a quantum state one may resort to full quantum 
state tomography [sj that, however, becomes impracti- 
cal in higher dimensions and may be affected by large 
uncertainty d, Other methods, requiring a reduced 
number of observables, are based on visibility measure- 
ments Bell' tests [H, [l^, entanglement witnesses 
or are related to Schmidt number [H-Illl. Mamr 
of them has been implemented experimentall y |23 -|27|. 
also in the presence of decoherence effects [28l. l29l| . 

As a matter of fact, any quantitative measure of en- 
tanglement corresponds to a nonlinear function of the 
density operator and thus it cannot be associated to a 
quantum observable. As a consequence, ultimate bounds 
to the precision of entanglement measurements cannot 
be inferred from uncertainty relations. Any procedure 
aimed to evaluate the amount of entanglement of a quan- 
tum state is ultimately a parameter estimation problem. 



where the value of entanglement is indirectly inferred 
from the measurement of one or more proper observ- 
ables "so] . An optimization problem thus naturally arises 
when one looks for the ultimate bounds to precision, i.e. 
the smallest value of the entanglement that can be dis- 
criminated according to quantum mechanics, and tries 
to determine the optimal measurements achieving those 
bounds. This optimization problem may be properly ad- 
dressed in the framework of quantum estimation theory 
fSlHi^ , which provides analytical tools to find the opti- 
mal measurement and to derive ultimate bounds to the 
precision of entanglement estimation. In particular, be- 
ing entanglement an intrinsic property of quantum states, 
we adopt local quantum estimation theory and look for 
optimal estimators maximizing the Fisher information 

In this paper, we address experimental determination 
of entanglement for two-qubit optical systems and apply 
quantum estimation theory to derive optimal estimators 
and ultimate bound to precision. This technique has been 
successfully applied in ^SS*] to estimate the entanglement 
of a pair of polarization qubit with the ultimate preci- 
sion allowed by quantum mechanics. Here we refine and 



extend the results of 35| in two directions: On the one 



•Electronic address: 'm.genovese@ii 
^Electronic address: giordaQisL iTP 



It] 



t Electronic address: , matteo. paris@fisica.unimi.it| 



hand we present a set of experiments aimed at estimating 
the amount of entanglement of a larger class of families of 
two-qubit mixed photon states. On the other hand, we 
have performed full state reconstruction using two dif- 
ferent tomographic sets of projectors in order to show 
explicitly that the evaluation of entanglement from the 
knowledge of the reconstructed density matrix provides 
a less precise determination. In our scheme entangle- 
ment, is evaluated through visibility measurements and 
estimators are built by a suitable combination of coin- 
cidence counts with different settings. Those estimators 
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turn out to be optimal and to provide estimation with the 
ultimate precision in the limit of Poissonian distribution 
of coincidence counts. In addition, we demonstrate ex- 
perimentally that optimality is robust against deviation 
from the Poissonian behaviour. Our approach allows en- 
tanglement estimation at the quantum limit, and it is 
also useful to compare different noise models using only 
information extracted from experimental data. 

The paper is structured as follows. In the next Section 
we briefly review the basic notions of local estimation 
theory, whereas in Section IlIII we apply them to estima- 
tion of entanglement of states belonging to two relevant 
families of mixed states. Section IIVI describes in details 
the experimental apparatus used to demonstrate our the- 
oretical results, which are described in the Section |Vl A 
detailed discussion of the experimental results is given 
in Section whereas Section IVlII closes the paper with 
some concluding remarks. 



II. LOCAL QUANTUM ESTIMATION THEORY 

We now give the basis ingredients for the local estima- 
tion theory starting with the classical case. Suppose we 
have a set of parameters A = (Ai, • • • , A„) G A C R" la- 
belling different states of the physical system of interest. 
A statistical model of our system is a set of probability 
distributions S = {px{x)\X £ A) such that n is the sam- 
ple space of the random variable x. The fundamental 
question in estimation theory is how to optimally esti- 
mate the unknown true values of the parameters A given 
a sequence of outcomes of measurement on the system 
{a;i, • • • , Xm}- From an geometrical information perspec- 
tive, this problem was first treated by Fisher who intro- 
duced for the case = 1 the now called Fisher informa- 
tion metric F{X): 



dxpx{x) di \ogpx{x) dj logpxix) 



dx 



dipx{x)djpxix) 
Px{x) 



(1) 



where di = d\.. F{X) is a positive definite matrix that 
represents a metric on the parameter space A and whose 
information geometric content is given by the best resolu- 
tion with which one can distinguish neighbouring points 
in the parameter space. The Fisher information met- 
ric is additive, therefore for a sequence of independent 
and identically distributed measurements with outcomes 
{xi, • • • , Xm}, F^''[X) = MF{\). The next step in the es- 
timation theory requires the introduction of the concept 
of estimator; the latter is any algorithm or rule of infer- 
ence, which allows one to extract a value for the unknown 
parameters on the basis of the sole knowledge acquired 
via the measurement process, i.e. the sequence of out- 
comes {xi, • • • ,Xm}- We say that the random variable 
A : — 7> A is an unbiased estimator if E[X\ = A i.e., 
its expected value coincides with the true value of the 



parameter (s). The ultimate bound on the precision with 
which one can estimate the parameters A is given by the 
Cramer-Rao theorem , which can be stated in terms of 
the covariance matrix Cov[A]y = E[XiXj] — E[Xi]E[Xj] 
as: 



Cov[A] > j^F{xy 



(2) 



In particular, for a single parameter the inequality reads 



Var[A] > 



1 



MF{X) ' 



i.e. the variance of the estimator, and therefore the preci- 
sion of any estimation procedure, cannot be smaller than 
the inverse of the Fisher information times the number 
of repeated measurements. In the general case, the in- 
equality for the variance of each of the parameters, i.e. 

Ya.r[X,]>j^[F{X)-%, 

holds only at fixed values of the others parameters. 

The previous results can be extended to the quantum 
realm, also taking into account all the possible measure- 
ments that one can implement on the systems. The quan- 
tum statistical model is given by a set of density opera- 
tors depending on the parameters A: S = {px\X €z A}. 
A measurement corresponds to a Positive Operator Val- 
ued Measure (POVM), i.e. a set of positive opera- 
tors £ = {Ei} such that J^i-^i-^l — ^ ^^'^ such that 
Pxii) = Tr[£'jpA] is the probability of having the z-th out- 
come. The Fisher information matrix F£(X) in Eq. ([!]) 
for a specific measurement process £ can then be writ- 
ten in terms of the classical probabilities pxii)- What is 
now specific to the quantum estimation process is that 
the optimization over measurement processes £ may be 
carried out. The problem has been solved in terms of the 
inequality {A > B means that A — B is a positive matrix) 



F£iX)<HiX) 



(3) 



that states that the Fisher information of any measure- 
ment process is upper bounded by the Quantum Fisher 
information H{X) (QFI). The latter is an n x n positive 
definite real matrix which can be expressed in terms of 
a set of n positive, zero mean operators called symmet- 
ric logarithmic derivatives (SLD) Li, each satisfying the 
following partial differential equation 
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dipx = 7^{Lipx + p\Li) 



(4) 



In particular, if one expresses the density matrix in its 
spectral decomposition 

PA = IV'iXV'il , (5) 

i 

the SLD pertaining to the z-th parameter is 

_ Pn+Pm 
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where 

diPX = ^ d^Pn |'0n)(V'n | + (7) 
n 

+ ^Pn(\di1pn){'4>n\ + |V'n)(9jV«|) , 
n 

accounts for the dependence of both the eigenvalues and 
the eigenvectors on the set of parameters A. In terms of 
the LiS the elements of the QFI can be written as: 



PA 



LjLi 



(8) 



By using the spectral decomposition of px, the QFI can 
be expressed in terms of the partial derivatives of the 
eigenvalues and of the eigenvectors as: 



{diPn){djP„ 



Pn 



Pn H" Pm 



(9) 



III. 



n 



ESTIMATION OF ENTANGLEMENT FOR 
TWO-QUBIT SYSTEMS 



We now apply the formalism described in the previous 
Section to obtain explicitly the ultimate bound to preci- 
sion on the estimation of entanglement for two relevant 
statistical models, i.e. for two families of two-qubit states 
that will be used in the following. 



A. The decoherence model 

The first statistical model we are going to deal with 
corresponds to the set of the states described by the fol- 
lowing two-parameter family of density operators 



g = p\^M + {l-p)D, 



where 



IV) = V9|HH> + yT^|VV) 



(10) 



(11) 



represents a pure polarization two-photon state with 
horizontal H and vertical V polarization, and D = 
q |HH)(HH|-t-(l-9) I VV)(VV I describes a mixed contribu- 
tion coming from the decoherence of \^p), p £ [0, 1]. We 
will refer to this set as the decoherence model for \^p). 
For the state g, both the two non zero eigenvalues 



A± = (1 ± ^l-4(l-p2)g + 4(l 
and their respective eigenvectors 

v± = ^={~f±{p,q),0,0,g{p,q)} 



P'^W) 



N± 

f±{p,q) 
9{p,q) 




(12) 



/g'^ip, q) ± fl{p,q) 
l-2q± v/l-4(l- 



p2)g- 



4(1-^2)^2 



depend on the parameters p, q. The straightforward cal- 
culations of the partial derivatives in Eq. ^ show that 
both the eigenvalues and the eigenvectors contribute to 
the diagonal and off-diagonal terms of the QFI. However, 
the sum of the different contributions results in a simpli- 
fied expression, and the QFI 



H{p, q) = diag 



4(1 -g)g 
l-p2 ' 



1 



q-q 



(13) 



is diagonal. From this expression we see that the variance 
on any estimator q for the parameter q is independent 
on the mixing parameter p and is bounded, apart from 
the statistical scaling, by the inverse of corresponding 
element of the QFI matrix 



Var[g] > 



g(l-g) 
M 



The lower bound is maximal in correspondence of g = 
1/2, i.e. when the state jV') is maximally entangled. We 
are now interested in estimating the value of entangle- 
ment of the overall state g. To this aim we remind that 
the negativity of entanglement defined as 



11^?' 



1 



(14) 



is a good measure of entanglement for two qubit systems. 
In Eq. (1141) Ta denotes partial transposition with respect 
to system A, and ||...||i is the trace norm. Entanglement 
negativity for states belonging to the decoherence model 
is given by 



e = 2p^il-q)q. 



(15) 



In order to reexpress the QFI in terms of the nega- 
tivity we make the change of variable p — >■ p, q — >■ 
{p — y^p2 _ e)/2p; the QFI changes according to the Ja- 
cobian of the transformation and the lower bound to the 
covariance matrix of the estimators p, e now reads: 



Cov[p,e] >ff-i(p,e) 

pHl-p')e 
p{l 



(16) 



p')e 



P{l-P^)< 



1 



From this expression we see that the lower bound for the 
variance of any estimator e of the negativity of the state 
g is independent on p and is minimal in case of maximal 
entanglement 



Var[e] > ^(1 



a. 



(17) 



B. The Werner model 



A second statistical model of interest for our analysis 
corresponds to the set of states described by the following 
two-parameter family of density operator 



P 



1-p 



1 (g) 1 . 



(18) 
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The states of Eq. (TTSl) are obtained by depolarizing the 
pure entangled state We will refer to this family as 
the Werner model for As in the previous example 
upon varying the parameter p we may tune the purity of 
the state, whereas the amount of entanglement depends 
on both parameters. The eigenvalues of g' depends only 
on p, whereas the eigenvectors depends only on q. The 
QFI matrix is thus given by the diagonal form 



H{p, q) = diag 



p- 



l + (2-3p)p' q{l~q){l+p) 



(19) 



and the inverses of the diagonal elements correspond to 
the ultimate bounds to Var(p) and Ya.v{q) for any estima- 
tor of p and q, either at fixed value of the other param- 
eter or in a joint estimation procedure. Entanglement of 
Werner states may be evaluated in terms of negativity. 



e — max < 0, 



1 



p (l + 4^9(9 - 1)) -1 



(20) 



which implies that Werner states are entangled for 



[1 + 4^9(1-9)]"' <P<1. 



Upon inverting Eq. (|20|) for p or q we may parametrize 
the Werner states using {p, e) and evaluate the QFI ma- 
trix H{p, e), their inverses and, in turn, the corresponding 
bounds to the precision of entanglement estimation. The 
main result is that the ultimate bound to the variance, 
depend only very slightly on the other free parameter (q 
or p). In other words, estimation procedures performed 
at fixed value of p or q respectively show different preci- 
sion, but the differences are negligible in the whole range 
of variations of the parameters. We do not report here 
the analytic expression of the inverse QFI at fixed p or q, 
which is quite cumbersome. However, as it can be easily 
checked, we note that the bound on the variance on e 
that can be derived by the expression of H{p,t)^^ sim- 
ply coincides to first order with the bound in Eq. p7)) 
already evaluated for the decoherence model. We there- 
fore use in the following, also for the Werner model, the 
bound given in Eq. p7)) . It can be shown that, for the 
set of values of p that will be relevant for our experimen- 
tal analysis, this approximation is negligible with respect 
to all the other sources of uncertainty. 



IV. EXPERIMENTAL APPARATUS 

The family of entangled states, investigated in our 
work, is constituted by polarization entangled states of 
the field obtained by coherently superimposing two or- 
thogonally polarized type-I parametric downconversion 
emissions (PDC), as schematically depicted in Fig. [T] 
The linear horizontal polarization of an argon laser beam, 
at wavelength A = 351 nm filtered by dispersion prism 
and Glan-Thompson prism (GP), is rotated at angle (j) 
by using half-waveplate (WPO). It is fundamental for 



our application that only the laser line A = 351.1 nm is 
used. For this reason we have introduced in the setup a 
prism as wavelength selector for eliminating wavelengths 
other than A = 351.1 nm. In particular the closest one at 
A = 351.4 nm, which could realize an unwanted phase- 
matching condition in our PDC setup. Then, the laser 
beam is addressed to a pair of non-linear beta barium 
borate (BBO) crystals (? = 1 mm), having optical axis 
in orthogonal planes, where PDC process occurs, result- 
ing in creation of biphotons with orthogonal polariza- 
tion [1^, [s^] • Upon changing the polarization of the UV 
pump, we change the amount of PDC light, generated by 
each crystal. For example, PDC occurs only in crystal 
one if the polarization of the pump beam is horizontal, 
while for having a balanced PDC process in both crystals 
we have set the angle 4> at 45°, having diagonal polariza- 
tion of pump beam. 

In order to compensate phase shifts, due to ordinary 
and extraordinary path in the crystals, we tilt the quartz 
plates QP, introduced between the halfwave plate WPO 
and BBO crystals, at angle thus fixing the relative 
phase between biphoton components generated in first 
and second crystal. 




D1 IF WP1 BS 



FIG. 1: (Color online) Experimental setup to generate po- 
larization entangled two-photon states with variable entan- 
glement and to estimate its value with the ultimate precision 
allowed by quantum mechanics. A continuous wave Argon 
pump laser beam with wavelength A = 351.1 nm is filtered 
with a dispersion prism and then passes through a Glan- 
Thompson prism and a half-wave plate WPO that rotates the 
polarization by an angle cf>. PDC light is generated by two 
thin type-I BBO crystals (^ = 1 mm). After the crystals the 
pump is stopped by a filter (UVF), and the biphoton field 
is split on a nonpolarizing 50-50 beam splitter (BS). Then it 
passes through half- wave plates (WPl, WP2) and interference 
filters (IF), centered at the degeneracy 702 nm. Finally the 
biphotons are focused on commercial single photon detectors 
(Dl, D2). 



In order to maintain stable the phase-matching condi- 
tions, BBO crystals and QP are placed in a closed alu- 
minium box internally covered by polystyrene used as 
thermic insulator. The box is equipped with a controlled 
heating system with a standard feedback circuit. We 
have experimentally verified that the temperature stabi- 
lization system ensures appropriate control on the phase 
shift. After the box the pump is stopped by an ultra- 
violet filter (UVF), and the biphoton field is split on a 
non-polarizing 50-50 beam splitter (BS). With the post- 
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election performed by a coincidence count circuit (CC), 
we can refer to our state as an optical ququart [s^ , which 
is entangled in two variables: polarization and spatial 
mode. 

In ideal conditions the output state is described by the 
pure state 

= cos<^|HH) +sin0e**('^)|VV) (21) 

where 0/2 is rotation angle of pump halfwaveplate WPO 
and corresponds to phase shift between pair of hor- 
izontal photons created in the first crystal and pair of 
vertical photons from the second crystal. After passing 
the half-waveplates (WP1,WP2) in each spatial mode, 
the biphoton field is projected into a linear vertical po- 
larization state by means of Glan-Thompson polarizers. 
Phase plates WPl and WP2 are mounted on precision 
rotation stages with high resolution and fully motor con- 
trolled, that allow rotating the polarization of the beams 
in the course of measurement process. Spectral selec- 
tion is performed by interference filters (IF) with central 
wavelength A = 702 nm and FWHM = 3 nm. Short focal 
lenses collimate resulting biphoton field into single pho- 
ton avalanche detectors (Dl, D2). Electrical signal from 
detectors is used by coincidence count scheme (CC) with 
time window r = 1 ns. 

The measurements performed at the output are de- 
scribed as projection of state into factorized linearly- 
polarized two-photon state: 

n.(a, /3) = |a + (a + ® 1/3 + s'^){(3 + s'^\ (22) 

where x = {s + 2s'}, s, s' — 0, 1. 




Quartz tilting angle (p, degree. 

FIG. 2: Probability of coincidence counts while performing 
projection measurement no(-f , f ) on state \tp^^) having = 
J as function of quartz plates tilting angle tp. 

In Fig. [5] we show the dependence of the probability of 
the coincidence counts 

po(<p) = (V'|^|no(J, J)|V'|^), 



as function of quartz plates QP tilting angle (p. The max- 
imum of this curve corresponds to phase shift between 
photon pairs ^{lpm) = and the output state is the Bell 
maximally entangled state 

l*+) = l^f^M>«|HH) + |VV), 

while the minimum of that curve corresponds to the max- 
imally entangled state 

I'i'") = H^^pJ « |HH) - |VV) . 

In this work we have fixed the tilting angle of quartz 
plates to have zero phase shift, thus, the family of states 
in Eq. (|21|) reduce to the one of Eq. (fTTj) where q = 

COS^ ((/)). 

V. ENTANGLEMENT ESTIMATORS 

In order to estimate the entanglement content of the 
states produced by the experimental set up described in 
the previous Section, one has to choose an estimator e to 
extract the value of entanglement from the experimental 
data. We will compare three different approaches: two 
are based on full tomography of the polarization two- 
photon and one is based on implementing the optimal 
estimator able to saturate the ultimate bound derived 
via the QFI. 

Quantum state tomography is an experimental pro- 
cedure providing full density matrix reconstruction of 
a quantum system. This is realized by means of a set 
of measurements performed on an ensemble of identical 
quantum systems Q. For a quantum state belonging 
four-dimensional Hilbert space at least 16 linearly in- 
dependent measurements are needed to reconstruct full 
density matrix and, typically, each measurement cor- 
responds to a local projection of the input two-qubit 
state. To be able to perform this set of 16 linearly in- 
dependent measurements we added a quarter-waveplate 
in each measurement arm just before the half-waveplates 
(WP1,WP2). The first used tomographic protocol (J16) 
[3 3 . koj l involves projective measurements performed di- 
rectly on some components of the Stokes vector. In par- 
ticular, the measurement set corresponds to projection 
onto polarizations HH, HV, VV, VH, RH, RV, DV, DH, 
DR, DD, RD, HD, VD, VL, HL, RL, where H, V, R, 
L, D, denotes horizontal, vertical, right and left circular 
and 45° diagonal polarizations, respectively. Here, for 
example, the measurement setting HR means measur- 
ing horizontal polarization on the first qubit and right 
circular polarization on the second qubit. Another ap- 
proach [4l|, \^ involves local projection of each qubit 
symmetrically placed on Poincare sphere. Extension of 
this method to four-dimensional case (R16) allows ob- 
taining higher fidelity of the reconstructed states [i^ Q 
with respect to the previous one. Once the density ma- 
trix of the generated state has been reconstructed, the 
negativity of the state can be evaluated inserting the re- 
constructed matrix elements in Eq. (|14p . The precision 
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the tomographic estimation of entanglement is limited 
by the uncertainties on the matrix elements. The overall 
uncertainty on the estimated value of entanglement may 
be evaluated by error propagating. In the following, after 
describing the implementation of optimal measurement, 
we will compare its precision with that of tomographic 
estimation. 

We first start to briefly describe the estimator for the 
class of states defined by Eq. ((TO)) . As already described 
in (ssj , an optimal estimator of the entanglement can be 
found by noticing that the expressions of the probabilities 
Pa;(e; a, /3) = Tr[ginj:(a, /?)] obtained by the projection of 
the state g on measurement operators in Eq. (1221) with 
X = 0, 1, 2, 3, allows writing the following set of unbiased 
estimators 



e(a,/3) 



V{a,l3) - cos(2q!) cos(2/3) 
sin(2Q:) sin(2/3) ' 



where V{a, /S) ~ po — pi — P2 + P3 is the expected value 
of two-qubit quantum correlations (QC). Furthermore, 
the estimators corresponding to the measurement angles 
a, P — ±7r/4 are optimal, as can be seen by evaluating 
the Fisher information 



p,j:{e]a,l3)[delogpa:{e;a,l3)f 



which for the chosen angles gives F^{j,j) equal to 
QFI. Then we have to express these optimal estima- 
tors, e = y(±7r/4, ±7r/4), in terms of the coincidences 
counts, which are the results of the measurement process. 
This can be done by fixing for example a — P — ^Tr/A 
and then, for each measurement run j — 1, .., Af — 40, 
one records the vector kj = {fcoj, fci.j, fc2j, fca^}, where 
kxj = fca;j(— 7r/4, — 7r/4), is the number of coincidence 
counts for the projector n^; defined in Eq. (l22t as mea- 
sured by the coincidence circuit during a single time win- 
dow of 10 seconds, and whose expected distribution is 
given 

Px{e;a,P) = Tr[gn^{a, P)] . 

Finally, we have to derive the probabilities 
— 7r/4, — 7r/4) in the expression of V{a,l3) in 
terms of the relative frequencies k,j;j{a, f3)/Kj, where 
Kj = J2x total number of coincidences. For 

large values of Kj the coincidence rates kx,j{a, 13)/Kj 
converges to the probability Px{€;a,/3). Therefore, the 
optimal estimator can be written as desired in terms of 
the coincidences' vector: e = e(kj). 

A second statistical model, which is a possible candi- 
date to represent the output of our experiment, is the 
Werner model of Eq. (|18p . From the physical point 
of view it corresponds to incorporate in our our scheme 
a portion of "fake" coincidences that results from dark 
counts of SPADs and from the influence of the ambient 
unpolarized luminescence. Since this light is unpolarized, 
its density operator can be described by the identity in 
(ITSl) . The distribution of coincidences is given by 

p'x{e;a,P)^Tv[g'Ilxia,P)], 



and the unbiased estimators for the mixing parameter 
and the entanglement negativity of the state by 

p' = ^(0,0) 

e' = -i + iy(0,0) + y(-7r/4,-7r/4). (24) 

where V{0, 0) — V{a = 0, (3 = 0) has been defined above. 
The estimators may be then written in terms of the coin- 
cidence vectors k^ , which was previously defined and that 
is used for F(-7r/4, 7r/4), and rj = {roj, ri j, r2j, rsj}, 
which is used in an analogous way to define the proba- 
bilities in for 1/(0,0) and whose elements are defined as 
rx.j = rx,j{0,0) i.e., the number of coincidence counts 



for the projector Ux 



with a 



0,/3 



0; in this 



case the total number of coincidences is Rj = X^x^^-j- 
(23) The estimators can then be written as p' — p'i^j), and 



e' = e'(kj,rj 



VI. RESULTS 



We first observe that for e(kj) and finite KjS the uncer- 
tainty in the estimation of the entanglement are mostly 
due to fiuctuations Skx in the coincidence counts kxj 
around their average values (kx) = J2j kx,j/M . Thus, if 
we want to establish under which conditions on the fiuc- 
tuations 5kx the variance of the estimator e(kj) satisfies 
the required bound, we have to implement standard un- 
certainty propagation with the derivatives dx = d/dkx 
evaluated for kx = {kx)j and assuming independence 
among fiuctuations at different angles, we have 



Var(e) 



E \9.e\'Skl 

X 

■[{{ko) + {k:i))\5kl+5k 



{k2)f {5kl + 5kl 



(25) 



If we now assume that the counting processes have a 
Poissonian statistics, i.e. Sk'^ = Var{kx) — (fc^)^, then it 
is straightforward to prove that 



Var(e) = 



{KY 



(fco + fc3)(fcl 



k^) = _L (1 _ £-2) 



i.e. QC measurements allow for optimal estimation of 
entanglement with precision at the quantum limit. Since 
the inverse of QFI is given by = 1 — for a 

wide range of two-qubit families of states [sO] , the above 
calculations suggest that this is a general result. In par- 
ticular, following the discussion at the end of section ITTll 
the above result is true also for the Werner state. In 
other words, given a source emitting polarization two- 
qubit states with coincidence counting statistics satis- 
fying the Poissonian hypothesis, then the experimental 
setup of Fig. [1] allows for optimal estimation of entangle- 
ment at the quantum limit by means of a QC estimator. 
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We finally note that in order to test the Poissonian hy- 
pothesis in our experiment we evaluated the Fano factor, 

2 

which is defined as F — — , where tr^ is the variance and 
fir is the mean of a random process in some time win- 
dow T. For a Poissonian process Fano factor should be 
equal to unity. In our experiment we had slightly differ- 
ent values [33|, but the method still allows for optimal 
estimation, thus showing the robustness of optimal mea- 
surement against deviation from Poissonian behaviour. 



A. Almost pure states 



parameter. By making use of the relations in Eq. (p4)) 
one can apply the same arguments to the Werner case and 
derive an appropriate expression for et- Upon evaluating 
the corresponding sample means and variances we can 
therefore obtain the first result of our analysis. This is 
illustrated on Fig. [3] where we report the estimated value 
of entanglement as a function of the reference one assum- 
ing, for the description of the output signals, the families 
g (left plot) and g' (right plot) respectively. Here the 
uncertainty bars denote the 3<t confidence interval and 
from this plots it is apparent that the experimental data 
are compatible with both models. 



The experimental setup of Fig. [T] allows for the prepa- 
ration of quantum states with high value of purity, 
namely having mixing parameter p close to unity. In 
these conditions both family of states in Eq.s (ITO)) and 
([T8| described in section HIT] are expected to give a reli- 
able estimation of entanglement. In order to verify this 
assessment, in the first part of our experiment we have 
performed measurements with different values of initial 
entanglement corresponding to different values of q, i.e. 
of the angle </> determined by WPO. We first consider 
the decoherence model of Eq. ([TU| . This model can be 
considered as a description of the decoherence mecha- 
nisms occurring in the experimental setup due to fluc- 
tuations of the relative phase between the two polariza- 
tion components, which results in fluctuation of phase 
shift between biphoton created in two crystals. Our ex- 
perimental procedure is based on M — 40 repeated ac- 
quisitions of coincidence vector kj — {koj,kij,k2j,k3j}. 
We have randomized the composition of kj over the se- 
quence of measurements to avoid spurious correlations, 
and finally we have estimated entanglement as the sam- 
ple mean (e) = £{kj)/M . The corresponding un- 
certainty has been evaluated by the sample variance 
Var{i)^j:Mk,)~{e)]yiM-l).^ 

In order to verify the compatibility of data with the 
decoherence model of Eq. pO)) we need to estimate the 
negativity with a second procedure, namely we make use 
of the estimation of the parameter p, quantifying the 
amount of mixing introduced by decoherence processes. 
We therefore define an unbiased estimator p b y first re- 
versing formula of the negativity i.e., p = ^e/y/q{l — q). 
We then note that the values of q and 1 — g in this 
model are given by the probabilities relative to the pro- 
jective measurements Ho (0,0) and 113(0,0) respectively, 
that can be expressed in terms of the elements roj and 
r3j . The estimator for p then reads 

Pir„k,)^\eik,)^^, 

where again Rj ^J2x ''^j - Rewriting the negativity de- 
fined in Eq. (1141) in terms of the pump polarization angle 
(j) we obtain e = psin2(/). Thus the reference value et of 
the negativity is then inferred as et = (p) sin 2(f>^ i.e. using 
the knowledge of and the estimation (p) of the mixing 



0.9 
0.7 
0.5 
0.3 



0.1 0.3 0.5 0.7 0.9 




FIG. 3: (Color online) Estimated value of entanglement as a 
function of the reference one assuming, for the description of 
the output signals, the families q (left plot) and q' (right plot). 
The uncertainty bars stays for the 3a confidence interval. 

Notice that the reference value is built, on the basis of 
a given model, in part with informations coming from the 
experimental settings (the tuning of the angle (j)) and in 
part from the results of suitably chosen coincidence mea- 
surements. On the other hand, the estimated value of 
entanglement is obtained solely with experimental quan- 
tities. In principle, we are not expecting the reference 
value to be more precise that the estimated one. The 
idea here is to use two different estimates of the same 
quantity (entanglement) obtained in two different and in- 
dependent ways in order to to discriminate and validate 
the different statistical models. Following our analysis, 
a given model is not suitable for the description of our 
system if the two different estimates that can be derived 
by that model, together with the resulting errors, are not 
compatible. 

It is interesting to compare these results, in particu- 
lar the ones which refer to the decoherence model (left 
plot in Fig. 121), with those obtained for a different set 
of measurements data presented in jssj . In that case a 
less precise control of the temperature of the PDC gen- 
eration system made more relevant the fluctuation of the 
phase and thus the state more mixed. Therefore, in that 
case, a self-consistent statistical analysis of the acquired 
data allowed discriminating between the two statistical 
models identifying the decoherence model of Eq. (fTO)) 
as the correct one for the experimental set up used in 
[35} . In the present case, which includes that the already 
mentioned control in temperature, the states obtained 
are nearly pure and thus one cannot expect the different 
characterization of noise to be relevant. Furthermore, to 
experimentally obtain more pure state one should reduce 
the collection angle of PDC emission. This obviously re- 
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duces the rate of coincidence counts, thus inducing an 
increase of the variance of both the estimators, for nega- 
tivity and purity parameter respectively. 



<e> <e> 




FIG. 4: (Color Online) Estimation of entanglement at the 
quantum limit. The plot shows the estimated value of entan- 
glement (e) according to the decoherence (left) and Werner 
(right) models as a function of the reference one et- The un- 
certainty bars on (e) denotes the quantity •\/Var(e) x {K), 
i.e. the square root of the sample variance multiplied by the 
average number of total coincidences {K) . The gray area cor- 
responds to values within the inverse of the quantum Fisher 

— 1/2 

information et ±iife J . Uncertainty bars on the abscissae are 
due to fluctuations in the estimation of the mixing parameter. 

We now pass to evaluate the optimality of our esti- 
mation procedure. In Fig. 2] we show, for the decoher- 
ence (left) and Werner (right) model, the estimated value 
of entanglement as a function of the reference one ob- 
tained for different values q = 0.97,0.93,0.88,0.78,0.5 
(i.e. = 10°, 15°, 20°, 28°, 45°). Note that the corre- 
sponding estimated mixing parameter in both model is 
larger than 0.97 for all points. The uncertainty bars on 
(e) denotes the quantity •\/Var(e) x (X), i.e. the square 
root of the sample variance multiplied by the average 
number of total coincidences {K) . This is in order to al- 
low a direct comparison with the Cramer-Rao bound in 
term of the inverse of the Fisher information (the gray 
area). Uncertainty bars on the abscissae correspond to 
fluctuations Set in the determination of et, due to fluc- 
tuations in the estimation of the mixing parameter with 
the procedure outlined above. The plot shows that our 
procedure allows estimating the entanglement with a pre- 
cision at the quantum limit for any value of q. From the 
figure it is also apparent that, due to the high purity 
achieved with the experimental set up that includes the 
active temperature control, and, therefore, due to the ir- 
relevance of the decoherence introduced, both the models 
give optimal estimation. Notice that this conclusion is 
robust against the fact that the statistics is not exactly 
Poissonian. 



B. Comparison with tomographic estimation 

We compared our results with estimation of entangle- 
ment from density matrix elements obtained exploiting 
two different procedures of quantum state tomography 
[39I [ioj . We found that the reconstructed density matri- 
ces are, for both tomography protocols, statistically com- 



patible within with both the two models of Eqs. and 
([T5|. As an example we present in Fig. [5] real and imagi- 
nary part of reconstructed density matrices of maximally 
entangled state corresponding to q = i (i.e., </) = 45°). 




FIG. 5: (Color online) Real (left) and Imaginary (right) part 
of the tomographically reconstructed density matrix for the 
maximally entangled state with J16 (top) and R16 (bottom) 
protocols. All the real elements, except the four dominant, 
and the imaginary ones are compatible with zero within the 
estimated tomographic uncertainties (not shown in the fig- 
ure). 



In fact, the tomographic procedure also allowed us to 
estimate entanglement and the corresponding variance. 
In order to have a fair comparison of the uncertainties ob- 
tained with different methods we have set measurement 
time for the tomographic reconstruction such to have thr 
total number of registered coincidences counts equal to 
M{K), i.e. the total number of coincidence in the op- 
timal measurement. The values of negativity calculated 
directly using the reconstructed density matrices and its 
variance (obtained by error propagation) for the maxi- 
mally entangled state are presented in Table U together 
with the determination obtained from the optimal mea- 
surement maximizing the QFI. All three negativity values 
overlap in their uncertainty intervals; the three methods 
are therefore coherent. Furthermore, it is evident from 
the presented results that the optimal method devised 
in this paper allows, at fixed sample size, for a sensitive 
reduction of the uncertainty in entanglement estimation. 
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Method 


e 




5e 


Optimal 


0.972 


± 


0.011 


Tomography: J16 


0.984 


± 


0.048 


Tomography: R16 


0.957 


± 


0.046 



TABLE I; Estimated value of entanglement with different 
methods. The uncertainty 5e is calculated usinv Eq. (|25|) 
for the optimal method and with error propagation for tomo- 
graphic estimation. 



C. Statistical mixtures 

In order to check our method in different working 
regimes we applied the estimation procedure to a set 
of mixed states obtained in a controlled way, i.e. by 
adding some portion of unentangled light to pure entan- 
gled state. As we have described in the previous section, 
our experimental set up allows us to obtain states with 
an extremely high purity. In the following we thus as- 
sume that the output state of our apparatus is the pure 
state as in Eq. ((TT|) . Then, if one is able, for example, to 
mix in a controlled way the components \HH){HH\ and 
|V^F)(V^F| to the maximally entangled states one obtains 
the states 



g = p\^m + {l-p)D 



(26) 



D = -{\HH){HH\ 



\VV){VV\)) , 



which correspond to the model ([TU|) with an adjustable 
mixing parameter. In practice, in order to tune the value 
of the mixing parameter p we have measured coincidence 
counts for states |HH)(HH| and |VV)(VV| for different time 
intervals. The sample of coincidence counts is then added 
to experimental data obtained for the maximally entan- 
gled pure state and then analyzed as in the previous sec- 
tion. 













0.1 0. 
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5 0.7 0.9 




FIG. 6: (Color Online) Estimation of entanglement at the 
quantum limit. The plot shows the estimated value of en- 
tanglement (e) as a function of the reference one tt- In the 
left panel we report estimated entanglement for mixed states 
generated according to the decoherence model (|26|) . In the 
right panel we report estimated entanglemed for mixed states 
generated according to the Werner model (|27|l . The points 
correspond to different portions of incoherent addition from 
both crystals. 



for the originally maximally entangled state (g = ^) 
and for states prepared with mixing parameter p = 
99.5%, 83%, 74%, 50%, 33%. 

A similar analysis may carried out for the Werner 
model. In this case, in order to tune the value of the 
mixing parameter p one should supplement the coinci- 
dences vectors kj and rj with values coming from un- 
polarized light. This can be achieved by measuring co- 
incidence counts for |HH)(HH|, |HV)(HV|, |VH)(VH| and 
|VV)(VV| for different time intervals. The measured val- 
ues are then added to the previously measured values for 
pure maximally entangled state. In this way, one can get 
data corresponding to 



(27) 



which correspond to a Werner state with tunable depo- 
larizing parameter. After performing measurement and 
analysis set described in previous section we can estimate 
entanglement and mixing parameter value in this family 
of states. In the right panel Fig. [5] we show the esti- 
mated value of entanglement as a function of the actual 
one for the originally maximally entangled state and mix- 
ture parameter p = 99.5%, 76%, 62%, 52%, 45%. As one 
can evince from the presented figure our method provides 
optimal entanglement estimation also for mixed states. 



VII. CONCLUSIONS 

In this paper we have addressed in detail the esti- 
mation of entanglement for pairs of polarization qubits. 
Our scheme is based on visibility measurements of quan- 
tum correlations and allows optimally estimating entan- 
glement of families of two-photon polarization entangled 
states without the need of performing full tomography. 
Our procedure is self-consistent and allows estimating the 
amount of entanglement with the ultimate precision im- 
posed by quantum mechanics. Although optimal estima- 
tion of entanglement does not require the full tomography 
of the states we have also performed state reconstruction 
using two different sets of projectors and explicitly shown 
that they provide a less precise determination of entan- 
glement. 

The technique has been demonstrated for nearly pure 
states as well as for controlled mixtures in order to con- 
firm its reliability in any working regime. With a suitable 
choice of correlation measurements it may be extended 
to a generic class of two-photon entangled states. The 
statistical reliability of our method suggests a wider use 
in precise monitoring of external parameters assisted by 
entanglement. 
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